The liquid to vapor phase transition in excited nuclei 
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For many years it has been speculated that excited nuclei 
would undergo a liquid to vapor phase transition. For even 
longer, it has been known that clusterization in a vapor carries 
direct information on the liquid- vapor equilibrium according 
to Fisher's droplet model. Now the thermal component of the 
8 GeV/c TT+^^^Au multifragmentation data of the ISiS Col- 
laboration is shown to follow the scaling predicted by Fisher's 
model, thus providing the strongest evidence yet of the liquid 
to vapor phase transition. 



Nuclear multifragmentation, the break up of a nuclear 
system into several intermediate sized pieces, has been 
frequently discussed in terms of equilibrium statistical 
mechanics, and its possible association with a phase tran- 
sition However, ample uncertainty remains regard- 
ing its nature, in particular whether multifragmentation 
is a phase transition and if so whether it is associated 
with the liquid to vapor phase transition. 

This paper will show that: 1) high quality experimen- 
tal data contain the unequivocal signature of a liquid to 
vapor phase transition through their strict adherence to 
Fisher's droplet model; 2) the two phase coexistence line 
is observed over a large energy/temperature interval ex- 
tending up to and including the critical point; 3) several 
critical exponents, as well as the critical temperature, the 
surface energy coefficient and the compressibility factor 
can be directly extracted; 4) the nuclear phase diagram 
can be constructed with the available data; 5) the nuclear 
liquid at break-up is a slightly super-saturated vapor, for 
which the pressure and density can be determined as a 
function of the temperature. 

In past attempts to investigate the relationship be- 
tween nuclear multifragmentation and a liquid to vapor 
phase transition, critical exponents have been determined 
l0'iF0 ' caloric curves have been examined |^ , and the ob- 
servation of negative heat capacities have been reported 



Q . Other studies of multifragmentation data have shown 
two general, empirical properties of the fragment multi- 
plicities called reducibility and thermal scaling P, p^ , |rT| . 

Reducibility refers to the observation that for each 
energy bin the fragment multiplicities are distributed 
according to a binomial or Poissonian law. As such, 
their multiplicity distributions can be reduced to a one- 
fragment production probability according to a binomial 
or Poissonian distribution. 

Thermal scaling refers to the feature that the aver- 
age fragment yield {N) behaves with temperature T as 
a Boltzmann factor: {N) oc exp(— B/T). Thus a plot of 
In (iV) vs. 1/T, an Arrhenius plot, should be linear. The 
slope B in such a plot is the one-fragment production 
"barrier" . 

Both the features of reducibility and thermal scaling 
are inherent to any statistical model, in particular to the 
clusterization of droplets from a vapor as described by 
Fisher's droplet model [|l2|. Thus it is interesting to see 
if a system portraying reducibility and thermal scaling 
portrays also the scaling of Fisher's model Fisher's 
droplet model describes the aggregation of molecules into 
clusters in a vapor, thus accounting for its non-ideality. 
The abundance of a cluster of size A is given by: 
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where ua is the number of droplets (or fragments) of 
mass A, normalized to the size of the system Aq; qq is 
a normalization constant depending only on the value 
of T [Q; T is the topological critical exponent; Afi = 
fi — fii, and /i and fii are the actual and liquid chemical 
potentials respectively; c^eA'^ is the surface free energy 
of a droplet of size A; cq is the zero temperature surface 
energy coefficient; a is the critical exponent related to the 
ratio of the dimensionality of the surface to that of the 
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volume; and e = (Tc — T)/Tc is the control parameter, a 
measure of the distance from the critical point, Tj,. 

Equation (|^) shows clearly that thermal scaling is con- 
tained in this description. In particular, at coexistence 
(A/i = 0), Eq. (0) leads to a "barrier" B = cqA" . 
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FIG. 1. (a) The ISiS fragment yield surface: natural log 

of the fragment yield vs. fragment mass and inverse tem- 
perature, (b) Arrhenius plots for representative charges, (c) 
Fragment mass yields for various values of E* . Solid curves 
are from a fit to Fisher's droplet model; dashed curves show 
the extrapolation of the fit results. See Fig. ^ for symbol def- 
inition. Error bars (statistical) are smaller than the size of 
the points. 

Recently, gold multifragmentation data from the ISiS 
Collaboration was shown to exhibit reducibility and ther- 
mal scaling in the fragment production probabilities 
[^,|6|. Since this behavior is inherent to Fisher's model, 
it is interesting to determine if Fisher's model describes 
the ISiS data set. In order to find if this is the case, the 
ISiS charge yields from 8 GeV/c ti^+ ^^^Au fragmenta- 
tion data (see Fig ^) were fit to Eq. (H). 

The mass of a fragment A (estimated by multiplying 
the measured fragment charge Z by the A-to-Z ratio of 
the fragmenting system) was used as the cluster size in 
Eq. (j^). The total number of fragments of a given size 
Na was normalized to the size of the fragmenting system 
Aq thus UA = Na/Aq. Here ^/E* was used in lieu of T, 
assuming the system behaves as a degenerate Fermi gas. 



The parameters of Fisher's model, e.g. t, a, and E* 
(in lieu of Tc) were used as fit parameters. The distance 
from equilibrium An was parameterized by a polynomial 
of degree four in E* and the coefficients of that polyno- 
mial were used as fit parameters. The surface energy co- 
efficient Co was parameterized by a polynomial of degree 
one and the coefficients of that polynomial were used as 
fit parameters. The level density parameter was absorbed 
into the fit parameters for A/x and Cq. 

While analyses similar to this one have been performed 
on multifragmentation data in the past p7|~pj|| , those ef- 
forts dealt with inclusive data sets; data from every exci- 
tation energy were examined as a whole and the results 
were presented as a function of incident beam energy. 
This work makes use of the high statistics, exclusive data 
set of the ISiS Collaboration and bins the events in terms 
of reconstructed excitation energy [Q. In addition, ex- 
plicit use of Fisher's expressions for the bulk and surface 
energies allows A/i and cq to be determined by the data. 
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FIG. 2. The scaled fragment yield distribution versus the 
scaled temperature. ISiS data (upper) from a wide range 
in charge collapse to a single curve over more than five or- 
ders of magnitude. Also shown are similarly scaled data 
from a three-dimensional Ising model calculation (lower); 
the quantity (nyi/goA~^)/100 is plotted against the quantity 
3.5(A''e/r) for clusters of size A = 5, 10, ... , 95, 100. 

Data for E* < E* and for 5 < Z < 16 were included in 
the fit to Eq. (0) and the parameters were allowed to float 
to minimize xZ- Fisher's parameterization of the surface 
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energy CqE is invalid for T > Tc, thus excitation energies 
greater than E* were not considered. Fisher's model ex- 
presses the mass/energy of a fragment in terms of bulk 
and surface energies. This approximation is known to fail 
for the lightest of nuclei where structure details (shell ef- 
fects) dominate the mass. For this reason and the fact 
that for the lightest fragments equilibrium and nonequi- 
librium production cannot be clearly differentiated, frag- 
ments with Z < 5 were not considered in the fit. Frag- 
ments with Z > 16 were not elementally resolved 
and were also excluded. 
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FIG. 3. The behavior of A/.1, and of the surface energy 

coefficient as a function of excitation energy up to = 3.7 

MeV/nucleon. See text for details. 

The behavior of the data for the {nA,A,T) surface is 
reproduced over a wide range in excitation energy and 
fragment charge. This is also shown in Fig. |^. In this 
figure the results of the analysis are presented in terms 
of Arrhenius plots (Fig. Qb), or of the fragment yield 
distribution plots (Fig. |l|c). A powerful method to ob- 
serve the results of this analysis directly is to scale the 
data according to Eq. (^ using the parameters result- 
ing from the fitting procedure. Figure ^ shows such a 
result. The fragment mass yield distribution is scaled 
by the Fisher's power law pre-factor and the bulk term: 
nA/qoA~'^ exp{AfiA/T). This quantity is then plotted 
against the temperature scaled by Fisher's parameter- 
ization of the surface energy: A'^e/Ve*. The scaled 
data collapse over five orders of magnitude onto a sin- 
gle curve, which is precisely the behavior predicted by 
Fisher's droplet model This curve is equivalent to a 
liquid-vapor coexistence line, as will be shown below, and 
provides the best, most direct evidence yet for a liquid 
to vapor phase transition in excited nuclei. 

To illustrate the generality of this type of scaling. Fig. ^ 
also shows the scaled cluster distributions from a three- 
dimensional Ising model calculation This system is 
know to model liquid-vapor coexistence up to the criti- 
cal temperature. The perfect scaling of the cluster yields 
according to Eq. (|^) demonstrates liquid- vapor-like coex- 
istence up to Tc- 

The value of r extracted from this analysis, r = 




T/T 

c 

FIG. 4. (a) The reduced pressure-temperature plot and (b) 
the reduced density-temperature plot. Filled symbols show 
the path of the ISiS data, open symbols show the coexistence 
curve of charged finite nuclear matter. 



2.00 ± 0.01, is in the range predicted by Fisher's model. 
Similarly, a = 0.64 ± 0.05 is close to the value expected 
for a three dimensional system, ~ 2/3. The values of 
A/i and cq and their energy dependence are shown in 
Fig. H The positive value of Afj, for E* < E* indicates 
that the system behaves as a slightly super-saturated va- 
por. Values for A/x and cq can be determined by making 
a choice of the level density parameter. For a = Aq/S, 
Co « 10.0±0.3 MeV while a = Ao/13 gives cq ~ 12.8±0.4 
MeV; both values are in the range of the value of the sur- 
face energy coefficient of the semi-empirical mass formula 
(^ 17 MeV), although the linear extrapolation all the 
way down to T = implied by Eq. (^) is questionable. 
The values of the critical exponents determined here are 
in agreement with those determined for the EOS multi- 
fragmentation data H] and the value of the excitation en- 
ergy at the critical point E* = 3.70 ± 0.05 MeV/nucleon 
is in the neighborhood of the value observed in the EOS 
analysis {E* = 4.75 MeV/nucleon) §,|J2|]. The EOS 
analysis relied on the assumption that A/i w 0. No such 
assumption was made in this work. In fact the present 
effort tests that assumption and finds it to be approxi- 
mately valid below the critical point since the observed 
A/i values are at most ~ 1 MeV (assuming a — Aq/IS). 
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Using the value of E* determined here and a level den- 
sity parameter of the critical temperature can be 
estimated to be on the order of 6.9 MeV, which is compa- 
rable to theoretical estimates for small nuclear systems 
1^,^. In most theoretical calculations the Coulomb 
force and small size of the system drastically reduces the 
value of Tc- It is also well known that both Tc and pc 
scale as a function of system size in many systems p5| , p6[ . 

The coexistence line and actual position of the 
fragmenting system in the pressure-density-temperature 
(P, p, T) diagram can be determined from this analysis. 

Fisher's theory assumes that the non-ideal vapor can 
be approximated by an ideal gas of clusters. Accordingly, 
the quantity ha is proportional to the partial pressure of 
a fragment of mass A and the total pressure due to all of 
the fragments is the sum of their partial pressures: 
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The reduced pressure is then given by: 

This is a transformation of the information given in Fig. ^ 
onto a more familiar frame of reference. 

When values of ua corresponding to A/i = are used 
in these equations, the coexistence curve is obtained. In 
other words, this analysis provides simultaneously the 
coexistence line and the actual line. Both arc shown in 
Fig §1. 

The system's density can be found via 
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As before, both the actual curve and the coexistence 
curve can be determined. Both are shown in Fig. Bb 
together with the coexistence curve of Guggenheim [E7|. 

The coexistence curve from the nuclear data does not 
agree closely with the Guggenheim plot over the full 
range of T/Tc. This is not surprising given the com- 
plexity of the nuclear fluid compared to the simple fluids 
analyzed by Guggenheim. However, near Tc, the univer- 
sal behavior of the coexistence curve is recovered. 

The compressibility factor Pc/TcPc is calculated by di- 
viding the denominator of Eq. (|^) by that of Eq. (^) and 
found to be 0.24 ±0.06. This value agrees with the value 
of the compressibility factor of many fluids ||2^ . 

In conclusion this paper has shown that the data of the 
ISiS Collaboration contain the unequivocal signature of 



a liquid to vapor phase transition via their strict adher- 
ence to Fisher's droplet model. Through Fisher's scal- 
ing of the fragment yield distribution (Fig. ||) the two 
phase coexistence line has been determined over a large 
energy /temperature interval extending up to and includ- 
ing the critical point. The critical exponents r and a as 
well as the approximate value of Tc, the surface energy 
coefficient cq and the compressibility factor have been 
extracted and agree with accepted values. A portion of 
the nuclear phase diagram has been constructed with the 
available data. Finally, the nuclear liquid at break-up 
was observed to be a slightly super-saturated vapor for 
which P/Pc and p/pc were determined as a function of 
T/Tc. 
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